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We calculate nuclear matrix elements (NME) of neutrinoless double beta decay {Oi/PjS) in four 
different candidate nuclei C^Ge, *^Se, ^"^"Mo, ^'^''Tc) within the quasiparticle random phase approx- 
imation (QRPA) and its uncertainties. We assume (up to) four coexisting mechanisms for 01/2/3 
decay, mediated by light Majorana neutrino exchange (v), heavy Majorana neutrino exchange (A*'), 
, i?-parity breaking supersymmetry (J?), and squark- neutrino (q), interfering either constructively or 

04 ■ destructively with each other. We find that, unfortunately, current NME uncertainties appear to 

prevent a robust determination of the relative contribution of each mechanism to the decay am- 
^ , plitude, even assuming accurate measurements of decay lifetimes. The near-degeneracy of Qv^p 

mechanisms is analyzed with simple algebraic techniques, which do not involve assumptions about 
the statistical distribution of errors. We discuss implications of such degeneracy on prospective 
searches for absolute neutrino masses. 
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<~i ' The discovery of neutrino masses and mixings 1 has given new impetus to the search for lepton number violation 
(LNV) and, in particular, for the process of neutrinoless double beta decay {0vl3j3) 0, 

; [Z, A) ^ {Z + 2, A)+2e- , (1) 

' in a variety of experiments using different [Z, A) nuclei [3]. The process can be mediated by light Majorana neutrinos, 
] as well as by alternative (and possibly co-existing) mechanisms invoking new particles and interactions beyond the 
(N ' Standard Model 

■ In the case of a single LNV mechanism (labelled by an index j), the Qvfifi decay half-life Ti in a given nucleus 
<^ , i — [Z, A) reads 

T-'^Gi\Mlri,\' , (2) 

> : ■ 

•'-j where Gl is a calculable phase-space factor, Mf is the Qu^P nuclear matrix element (NME), and r]j is a dimensionless 
/\ ■ LNV parameter, characteristic of the particle physics model. In the case of coexisting mechanisms with identical 
I phase space {G{ = Gi), the above expression is generalized as [51 
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T-' = GAY,M^rj, . (3) 



In the above equation, the rjj parameters may take either sign (and, allowing CP violation, even complex phases), 
leading to constructive or destructive interference in the decay amplitude. 

In the absence of independent constraints on the LNV parameters, the discrimination of different mechanisms 
requires that either the r]j's are factorized out, or that they are determined through the data themselves. For instance, 
assuming one mechanisms at a time, one may analyze 77-,- -independent ratios of half-life data {Ti} from different nuclei 
(typically with a ^^Ge datum at denominator) possibly including joint {Ti, Tk) probability distributions fl0| . 

Alternatively, one may assume (at least) as many {Ti} data as the number of considered mechanisms, and then solve 
Eq. ^ in terms of the r]j 's ^ . 

The studies in suggest that, at least in principle, a few (possibly coexisting Q) mechanisms may be dis- 

criminated with accurate half-life data, provided that the associated NME uncertainties in different nuclei are "small 
enough." Attempts to quantify this statement in prospective scenarios have been made through MonteCarlo simu- 
lations with educated error guesses [9|] or through likelihood analyses ^Oj] including error covariances pjj]. However, 
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apart from the difEculty in estimating "la" theoretical errors, the approaches in [9|, |10| had to rely upon a somewhat 
inhomogeneous comparison of NME's taken from different publications or nuclear models; conversely, the NME's in 
were all calculated in the same model, but their variations were not fully included in the analysis, although it was 
noted that they could be relevant. In this work we adopt an approach where: (1) all NME's and their variations are 
estimated with one and the same method, the quasiparticle random phase approximation (QRPA [12]) and (2) the 
resulting discrimination of different mechanisms is analyzed with simple algebraic (rather than statistical) tools. 

In particular, the effect of theoretical uncertainties is evaluated herein by varying several QRPA inputs, which 
require extensive NME calculations. As a consequence, our analysis is restricted to a relatively small set of candidate 
nuclei and mechanisms. We consider four among the experimentally promising 3] candidate nuclei, labelled by the 
index i: 



i^^^Ge, 82se, ioomo, i30Te, (4) 

and four possible Oi^/3/3 decay processes, labelled by the index j, 

J^v,N,^,q, (5) 

corresponding to light Majorana neutrino exchange [isl ]. heavy Majorana neutrino exchange i?-parity breaking 

supersymmetry mechanism [T3-0] , and squark- neutrino mechanism (2ll423j | , respectively. We then study solutions to 
Eq. ^ in terms of prospective {T^} data. Unfortunately, we find that the reconstructed LNV parameters r^j of the 
four mechanisms to decays are destabilized once different, currently admissible QRPA variants are considered. 

The effective degeneracy of mechanisms may have significant implications on the interpretation of upcoming 
results, and on their interplay with other observations sensitive to absolute neutrino masses. 

Our work is structured as follows. In Section II we describe the LNV mechanisms considered for decay. In 

Section III we discuss relevant aspects of the nuclear structure calculations. In Section IV we analyze the problem of 
discriminating the mechanisms within QRPA uncertainties. In Section V we summarize our results and discuss their 
implications for absolute neutrino mass searches. 



II. LNV MECHANISMS 



In this section we describe relevant features of the four LNV mechanisms in Eq. ([5]). In each i-th nucleus [Eq. (|3])], 
they share the same phase space G^, estimated in [16] as 

G(^^Ge, *^2se, ^^Mo, ^^^Te) = (7.93, 35.2, 57.3, 55.4) x lO-^^y-^ . (6) 

Conventionally, the above phase space factors embed the fourth power of the axial coupling constant at the fixed value 
gA — 1.25, and the inverse square of the nuclear radius R — roA^/"^, with fixed tq = 1.1 fm. They are compensated 
by factors ((?yi/1.25)^ and 1/i? in the definition of the matrix elements Mf [23 |. 



A. Light Majorana neutrino exchange (u) 



In the usual case of exchange of three light Majorana neutrinos Vh (with masses m/^, Majorana phases (j)^, and 
mixings Ueh)i the LNV parameter takes the well-known form 

'^-^ me nie 

h=l 

where rrie is the electron mass, and m^^ is the "effective Majorana mass." In each nucleus, the nuclear matrix element 
M'^ consists of Fermi (F), Gamow-teller (GT) and tensor (T) contributions. 




with F, GT, and T operators defined as in 
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B. Heavy Majorana neutrino exchange (A'^) 

The neutrino mass spectrum may include, in principle, any number of heavy Majorana states Nk with masses 
M/j above the typical energy scale, i.e., greater than 0(10) GeV. We restrict our consideration to the case of 

right-handed Majorana fermions, which are singlets of the 811^(2) x U(l) group. Then, heavy neutrinos can mediate 
Oi^/3/3 decay with a LNV parameter ijn given by 

heavy 

W = E |t^e.|V*^^, (9) 
k ^ 

where nip is the proton mass, while the Uek are the mixing matrix elements associated to left-handed weak interactions. 
Also in this case, the nuclear matrix element is of the form 




with F, GT, and T operators discussed at length in [TH|- 



C. R-parity breaking mechanism {fCj 

In supersymmetric (SUSY) models one may define a multiplicative quantum number R — j2S+3B+L ^-y^rith S, 
B, and L being the spin, baryon, and lepton numbers), which takes the value +1 for ordinary particles, and —1 for 
their superpartners. i?-parity breaking allows (among others) superpotential terms of the form Wjji 3 X'^jf.LiQjD'^, 
where L and Q are lepton and quark doublet left-handed superfields, D'^ is a down quark singlet superfield, and A' is a 
trilinear coupling. These terms can trigger 01/(3/3 decay with short-range exchange of heavy superpartners (see below 
in this subsection) and with long-range exchange of heavy squarks plus light neutrinos (see the next subsection). 

Assuming a dominant, short-range exchange of a heavy gluino with mass rrig, the LNV parameter rj^ takes a 
simplified form, symmetric in the masses of the u and d squarks: 



6 Gj, m~g ml^ ■ m\ 



(11) 



where Gp is the Fermi constant and as the strong coupling constant. At hadron level, in the hypothesis of one-pion 
and two-pion exchange dominance, the matrix element contains both tensor and Gamow- Teller contributions, 

= (-^) ^ [c^''{M\'' + M(^^) + c2^(a4^ + M'^r^)] , (12) 
V 1 .25 ' 

with known prefactors c^^ and c^'^ . See [l^, [2^ for further details. 



D. Squark- neutrino mechanism (g). 

i?-parity breaking SUSY may trigger Qvl3f3 decay via long-range exchange of (not necessarily massive) neutrinos 
plus squarks [27l |. The corresponding LNV parameter is sensitive to down-squark ((i(fe) — d, s, b) masses and left-right 
mixings: 

fc \ <ii(fe) <i2(fc) 

At hadron level, in the hypothesis of pion-exchange dominance, the matrix element contains both tensor and Gamow- 
Teller contributions, with a resulting matrix element of the form 

= (i^)'^^'^St-M^) . (14) 



See 271 for further details. 
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E. Remarks on v, N, ft, q mechanisms 



The above four mechanisms are different in several respects: hght and heavy neutrino exchange NME's do not exactly 
scale as g\ (because of Fermi contributions), unlike SUSY mechanisms, which also have large tensor contributions; 
moreover, short-range mechanisms (such as heavy neutrino exchange) are particularly sensitive to the nucleon-nucleon 
potential. Thus, one might expect to find different NME isotopic patterns in different mechanisms. However, we shall 
see in Sec. IV that the resulting differences in NME patterns are comparable to current QRPA variations, thus 
preventing an effective discrimination of the mechanisms based on multi-isotope half-life data. 

The similarity of NME patterns in the different nuclei and mechanisms considered herein suggests that nucleon 
(rather than nuclear) physics dominates in all cases. For the v mechanism, the dominance of short-range nucleon- 
nucleon contributions has been shown in different approaches [25l[2S.i2i^]. It would be useful to perform similar studies 
for nonstandard mechanisms as well. 



In what follows, the Qufiji nuclear matrix elements Mf are evaluated (for i = ^^Ge, ®^Se, ^""Mo, ^'^°Te and j = 
V, N, fi, q) under the so-called self-consistent renormalized QRPA (SRQRPA) [30l - l32| . which takes into account the 
Pauli exclusion principle and conserves the mean particle number in correlated ground states. The following variants 
are introduced to account for intrinsic QRPA uncertainties. 

For each nucleus, two choices of single-particle basis are considered: intermediate and large. The intermediate-size 
model space (m.s.) has 12 levels (oscillator shells N = 2-4) for ^^Ge, ^^Se, 16 levels (oscillator shells N = 2-4 plus 
the f + h orbits from iV = 5) for ^""^Mo and 18 levels (oscillator shells iV = 3, 4 plus f + h+p orbits from N = 5) for 
i30rpg^ The large-size single particle space contains 21 levels (oscillator shells N — 0-5) for ^^Ge, ^^Se, ^°°Mo, and 23 
levels {N = 1-5 and i orbits from = 6) for ^'^'^Te. In comparison with previous studies [IJ] we omit the small-size 
model space, as it is insufficient to describe realistically the tensor matrix elements. 

The single particle energies are obtained by using a Coulomb-corrected Woods-Saxon potential. Two-body G- 
matrix elements are derived from the Argonne and the Charge Dependent Bonn (CD-Bonn) one-boson exchange 
potential within the Brueckner theory. The schematic pairing interactions are adjusted to fit the empirical pairing 
gaps [1^ . The particle-particle and particle- hole channels of the G-matrix interaction of the nuclear Hamiltonian H 
are renormalized by introducing the parameters gpp and gph, respectively. The calculations are been carried out for 
gph = 1.0. 

The particle-particle strength parameter gpp of the SRQRPA is fixed by the data on the two-neutrino double beta 
(2i//3/3) decays [2J|. In the calculation of the Oi^/3/3-decay NME's, we consider the two-nucleon short-range correlations 
derived from the Argonne or CD-Bonn nucleon-nucleon (NN) potentials as residual interactions fsS]. The measured 
2j//3/3-decay half-lives T2^, taken from 35], fix the matrix element |Mq^^| via Tj"^ = G2,yg\\meMQ^\'^ , where is 
the 5A-independent phase space factor for 2j//3/3 decay. The value of gpp is thus fixed for any given gA, which we take 
equal to either its bare value {gA = 1.25) or its quenched value {gA = 1.00) for all nuclei. For simplicity, we neglect 
the possibility of different gA quenching in different nuclei 36]. Finally, for each 2://3/3 constraint to gpp, we take both 
the T21, central value and its ±la deviations, leading to a triplet of gpp estimates (best fit ±lcr) in each case. 

Summarizing, our analysis is based on the calculation of 384 NME's, corresponding to four nuclei {i), four mech- 
anisms (j), two model space sizes (m.s.), two nucleon-nucleon potentials {NN), two axial-coupling values {gA), and 
three particle-particle strength values {gpp). 



However, wc have verified a posteriori that gpp variations within ±1(t induce NME variations much smaller than those 
induced by changes in model space, nucleon potential, or gA- The analysis discussed in the next Section turns out 
to be basically unaffected by such extra gpp variations. In what follows, we shall thus focus on the numerical results 
obtained at the central values for the gpp's. In particular, we restrict our analysis to eight 4x4 matrices M| of 
NME's, corresponding to eight QRPA variants, associated to changes in model space size, NN potential, and gA, 



III. NUCLEAR STRUCTURE CALCULATIONS 



384(NME) = 4(i) x 4(j) x 2(m.s.) x 2{NN) x 2{gA) x 3{gpp) . 



(15) 



8(M/) = 2(m.s.) X 2{NN) x 2{gA) ■ 



(16) 



Table U reports the corresponding M- estimates with three significant digits. 
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TABLE I: Matrices A'/f of NME's for i = (^^Ge, *^Se, ^°°Mo, ^^°Te) and j = {v, N, fi, q), for each of the 8 variants adopted 
within the QRPA. The 2x2x2 = 8 variants correspond to different choices for the model space size (intermediate or large), 
for the A'' A*' potential (Argonne or CD-Bonn), and for the qa value (1 or 1.25). 
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IV. DISCRIMINATION OF MECHANISMS: EXPECTATIONS AND RESULTS 

The system of equations in ([3]) involves, in general, four unknown complex parameters "qj. For the sake of simplicity, 
as in 0, we restrict ourselves to a CP-conserving scenario with real rjj. The relative sign of the ?7j's can, however, be 
positive or negative, leading to constructive or destructive interference among different mechanisms. Taking square 
roots in Eq. ([3]) leads to a further ± sign ambiguity for each index i [sj. We restrict our analysis to the case with all 
signs for each i, namely, to the linear system of equations 

4 

^Af/r;, (i = l,...,4) , (17) 
i=i 

where 

= (G.T,;)-^/' . (18) 
Relaxing the above simplifying assumptions would corroborate the results which follow. 
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A. Expectations 

For given multi-isotope data {T,}, solutions to the linear system in Eq. (jl7p in terms of {rjj} should identify, in 
principle, the underlying physics mechanism(s) [5j]. Testing different mechanisms, as well as their possible coexistence, 
becomes then a parameter estimation problem in 77j-space, provided that theoretical uncertainties (for M^) and 
experimental errors (for Ti) are accounted for. For instance, if Qv^jS decays were due only to light Majorana neutrino 
exchange, accurate solutions to Eq. ([T7|) should ideally cluster around a point (rj^, t/n, tj^, rjq) ~ (('7i/), 0, 0, 0). 
Conversely, nonzero values for two or more rji^ would signal the coexistence of two or more mechanisms in 
decay. This approach has been discussed in Q where, however, the NME values were taken from a single QRPA 
calculation with fixed input parameters, whose possible variations were not considered. 

As we show below, it turns out that the solutions {rij} are destabilized once QRPA variants are accounted for, 
even in the presence of accurate measurements (STi/Ti <C 1). In particular, the linear system of Eq. (jl7l) is shown 
to be ill-conditioned, with matrices Mf of nearly rank one, leading to an effective degeneracy of mechanisms. Before 
dealing with the M- algebra, it is useful to illustrate the source of such degeneracy in a graphical way. 

Figure 1 shows the contents of Tab. Has a scatter plot of NME's, for each couple of nuclei among ''^Ge, ®^Se, ^™Mo, 
and ^'^"Te. For each mechanism, the eight QRPA variants correspond to eight markers (circles for i/, triangles for N, 
squares for fi, and stars for q), solid and empty markers corresponding to = 1 and gA = 1.25, respectively. NME's 
for the v mechanism are multiplied by a factor of 100 in order to fit in the chosen logarithmic scales, which span the 
range [10^, 10'^]. Note that, on such scales, the effects of gpp variations mentioned in Sec. Ill would be unobservable. 
Diagonal dotted lines mark the band where the x and y coordinates are equal up to 20%. The source of degeneracy is 
then rather clear: all NME's cluster along the diagonal in each panel, with little spread (less than about 20%) in the 
orthogonal direction. In particular, the transverse spread induced by different mechanisms appears to be comparable 
to that induced by the QRPA variants. In other words, the columns of the M;^ matrices are roughly proportional 
to each other, and with comparable entries for different nuclei, within QRPA uncertainties. The linear system in 
Eq. P?)) is then expected to be significantly degenerate, with nearly "null determinant" not only for the full 4x4 
matrix , but also for all of its 3 x 3 and 2x2 submatrices. These expectations are confirmed by the following 
numerical arguments. 



B. Numerical results 



Stability and degeneracy issues in linear equation systems are best dealt with the Singular Value Decomposition 
(SVD) technique [U HI]- For a square matrix M as in Eq. (fT7|) the SVD is equivalent to diagonalizing the product 
matrix M • M-^, taking the square roots of the eigenvalues, and sorting them in descending order. The resulting 
set of "singular values" ai encodes relevant information on the stability of solutions [13, llM]- In fact, if all cr^'s 
have comparable magnitudes, the rjj solutions to Eq. (jl7p tend to be stable under small variations in either the 
Afj^ or fii entries. Conversely, if at least one cr^ is much smaller than the others, the system is ill-conditioned and 
generates unstable solutions. In particular, an ill-conditioned system, characterized by a large "condition number" 
K = Cmax/fmin, IS cxpccted to bccome unstable when the fractional variations 6Mf or S^i become comparable to 
the reciprocal condition number p = 1/k (stI . |38| . A slightly more technical point concerns the optimization of p. 
Although p is invariant under an overall scaling of matrix entries (Af| / • M|), it may change under separate 
scalings of rows and columns {M^ — )■ • M:^ ■ Sj), which transform Eq. (jl7p into an equivalent system, provided that 
rii —> ri i/si and pi = pi ■ si. The maximization of p by row/column rescalings, usually called "matrix equilibration" 
[37H39'| . provides more robust estimates of the fractional errors which may trigger unstable solutions. Summarizing, 
the SVD of an equilibrated matrix provides useful information about the effective matrix rank (equal to the number of 
"large" singular values) and about the size of "dangerous" fractional errors (equal to the reciprocal condition number 
P)- 

Table [III reports the singular values and (reciprocal) condition numbers, as obtained by applying equilibration and 
SVD routines 39] to the eight matrices in Tab. HI In each of the eight QRPA variants, only one singular value is 
of 0(1), the others being smaller by one to three orders of magnitude. The matrices are thus nearly of rank one, 
leading to basically undetermined and degenerate solutions. Instabilities of solutions to Eq. (ITT)) are expected to arise 
already for 6Ml as small as p ~ ©(lO"'^), implying that even three-digit approximations are potentially dangerous. 
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OulSiS NME scatter plot 
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FIG. 1: Graphical representation of the NME's in Tab.|T]as a scatter plot, for each couple of nuclei among ^^Ge, ^^Se, ^""Mo, 
and ^^°Te. Eight QRPA variants are shown for each of the four different Of (5/3 mechanisms, labelled as v (circles), A'^ (triangles), 
J4 (squares), and q (stars). Solid and empty markers refer to = 1 and qa = 1.25, respectively. In each panel, the {x,y) 
coordinates span the range [10^, 10^] in logarithmic scale. For the u mechanism, NME's are scaled (xlOO) to fit the range. 
Diagonal dotted lines mark the band where the x and y coordinates differ by up to 20%. The NME's appear to cluster along 
the diagonals, suggesting an effective degeneracy of the mechanisms. See the text for details. 
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TABLE II: Singular value decomposition of the (equilibrated) matrices Ml given in Tab.|Il For each of the eight QRPA variants, 
we report the four singular values Oi (in descending order), the matrix condition number k, and the reciprocal condition number 
p = The inequality a\ S> o"2,3,4 signals that each matrix Ml is ill-conditioned and nearly threefold degenerate. See the 

text for details. 
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Let us discuss a representative numerical check, by assuming a scenario with true NME's taken from the first 
matrix of Tab. HI and true LNV parameters (ryj/, rj^, ij^i, rjg) set to (6.18 x 10"^, 0, 0, 0), corresponding to a single 

QvjSp mechanism, namely, light Majorana neutrino exchange with mpp = 0.316 eV. Such choices for and rjj imply 
the following /i^ values (up to three significant digits): 

At(™Ge, *^2se, i°°Mo, ^^^Te) = (2.38, 2.22, 2.24, 2.03) x 10"^ , (19) 

and a ''^Ge half- life of 2.23 x 10^^ y, which purposely reproduces the claimed measurement in [10, Elj. With the above 
/ii's, we solve Eq. (|17p in the unknowns rjj for each of the eight matrices in Tab. HI including the first one (where we 
expect to recover the "true" rjj by construction). 

Table IIIII shows the results of our exercise, in terms of the relative contribution Wj of each j-th mechanism to the 
QvPP amplitude in ^^Ge (z = 1), 

W, = ^^ ^Y.W,=l. (20) 
Y^Mirj, 

In the first row of Tab. IIII[ one would expect to recover the true v mechanism, corresponding to Wi, = 1 and 
^ — 0. However, this is not the case: the three-digit approximation is already sufficient to generate spurious 
contributions (e.g., Wq ~ —0.16). In the second to eighth row in Tab. IIIII where QRPA variants entail relative 
variations SMf 3> 0(10-'^), disparate solutions emerge, with constructive or destructive interference of (sub)dominant 
mechanisms. For instance, in the eight row, three large coherent contributions from the ly, N, and ^ mechanism are 
almost cancelled by an even larger opposite contribution by the q mechanism. In Tab. IIIII the largest amplitude 
corresponds to Wi, in only two cases (1st and 3rd variant), while in others it corresponds to either Wn (5th variant) 
or Wq (remaining variants). 



TABLE III: Relative contribution Wj of the j-ih Qvj3P mechanism to the total ^^Ge decay amplitude (^ . Wj = 1), for eight 
QRPA variants, assuming Wj — 5ji as "true" case, and precisely known lifetimes Ti for the four nuclei (see the text for details). 



Variant 


W, 


Wn 


Wf, 


W^ 


1 


+1.05 


-0.01 


+0.12 


-0.16 


2 


+3.14 


+5.62 


-1.05 


-6.71 


3 


+1.26 


-0.24 


+0.94 


-0.96 


4 


+3.83 


+2.63 


+2.13 


-7.59 


5 


+1.00 


+3.07 


-2.83 


-0.24 


6 


+7.90 


+16.22 


+1.86 


-24.98 


7 


+1.97 


+1.09 


-0.27 


-1.79 


8 


+6.64 


+5.90 


+5.99 


-17.53 
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Note that the instabihty of the Wj's (in both magnitude and sign) arises from intrinsic properties of the NME 
matrices M| and would be present also for other choices of the "true" LNV parameters (e.g. by assuming another 
mechanism as dominant). Moreover, varying input data in Eq. (jlQI) within prospective errors, as well as removing the 
simplifying assumptions preceding Eq. ()17|) . would only generate further instabilities or multiplicities in the solutions. 

If the problem is reduced to just three or two mechanisms (instead of four), the instabilities, although less severe, 
remain harmful. In generic 3x3 and 2x2 minors of the Mil matrices, unstable solutions are expected to arise for 
fractional variations SMil of order = (Js/ai ^ few x 10^'^ and p2 ~ o-2/<^i ^ few x 10^^, respectively. Since current 
QRPA variants induce ^ O(10~^), a robust reconstruction of the "true" scenario cannot be achieved. 

As a numerical check, we have repeated the previous exercise in the subcase (''^Ge, ^■^'^Te) x {i^, q), which is among 
the "less degenerate" 2x2 ones. We find reciprocal condition numbers P2x2 — few %, as expected. Since QRPA 
variants entail ^ 0(10^^), solutions are unstable, and the reconstruction of true case {Wi, — 1) yields fake 

values for W^, ranging from a minimum of 0.23 (corresponding to dominant q mechanism) to a maximum of 1.5 
(corresponding to destructive interference of v and q mechanisms). If we further change either /i(^^Ge) or yLt(^^°Te) 
by ±10% (corresponding to ±20% experimental uncertainties in 0i//3/3 half-lives), we find reconstructed values of 
Wi, within a wider interval, Wi, S [—1.4, +2.9]. This range spans disparate solutions, including coexisting v and q 
mechanisms with constructive or destructive interference, as well as subcases close to single v or single q mechanism 
{Wi, ~ 1 or Wq ~ 1, respectively). Analogous results hold for other 2x2 minors. Note that, since the M/ matrices 
are nearly of rank one (cri ^ (J2.3a), the degeneracy applies to rectangular sub-matrices as well; thus, the situation 
would not significantly improved by over-constraining the system, e.g., by forcing two mechanisms to (approximately) 
fit more than two lifetime data among those considered. 

We stress that the observed degeneracy derives from the nearly threefold singularity of the matrices , and not from 
specific choices of prospective input data /i^ or of true LNV parameters rjj . The phenomenological indistinguishability 
of the {v, N, fi,, q) mechanisms in the (^^Ge, *^Se, ^°°Mo, ^■^'^Te) nuclei appears to be an unfortunate but general 
result, as far as the current QRPA framework and its estimated variants [leading to SMf ^ 0(10^^)] are assumed. 
The numerical analysis suggests a very ambitious accuracy goal for prospective QRPA estimates [i.e., a relative spread 
SM:^ at the few percent level], in order to clearly discriminate the relative contributions of two mechanisms by means 
of two nuclei, among those considered herein. 



C. Discussion 

The difficulty of separating the v mechanism from heavy Major ana neutrino exchange or from SUSY-mediated 
decays in multi- isotope Oi//3/3 data had been noted earlier in |8l-[l0l]. The results presented in this work corroborate 
such findings, with the advantage of being based on a homogeneous comparison of NME's, calculated in well-defined 
variants of one and the same nuclear model, for all the nuclei and mechanisms considered — which was not the case 
in [8l-[l0|. Therefore, the emerging degeneracy appears to be quite general and not accidental, in the context of 
the considered mechanisms: for any of them, the NME's turn out to have a similar size in the various nuclei, with 
uncertainties mostly absorbed by common scaling factors. The resulting numerical degeneracy in the determination of 
LEV parameters is supported and quantified by simple algebraic arguments, which do not assume specific distributions 
of theoretical errors. As such, the analysis method could be extended to other nuclei, mechanisms or nuclear structure 
models. 

Taken together, the results of [8l-[l0| and of this paper show, from different viewpoints, that 0i//3/3 mechanisms 
mediated by (light or heavy) neutrino exchange and by fi, supersymmetry may be phenomenologically degenerate 
in some nuclei of interest, within realistic uncertainties. Therefore, Oz^/3/3 signals or limits from such nuclei could be 
generically interpreted in terms of unconstrained linear combinations of coexisting {i^, N, fi, q) mechanisms, including 
subcases with one dominant (if not single) mechanism. Breaking the degeneracy would require significantly more 
accurate NME estimates than can be obtained in current, state-of-the-art QRPA. 

It should also be noted that we have considered only CP-conserving cases [i.e., real fii in Eq. (fT7|) ] and, furthermore, 
we have discarded the discrete ambiguities related to the two possible signs for each index i (/i^ — >■ ±/ii) [5], which 
would worsen the degeneracy problem by increasing the multiplicity of solutions. In CP-violating scenarios, with 
unknown relative phases (pi among different 0iy(3/3 decay amplitudes, the discrete ambiguities would become continuous 
(/ii — ^ e^'^^Hi), thus rendering the system in Eq. (jl7p undetermined, unless extra independent data were added to 
constrain the unknowns (f>i. In any case, removal of our simplifying assumptions about the CP properties of the fii 
parameters would make the degeneracy problem more (not less) difficult to solve. 

These negative findings do not exclude, of course, more favorable scenarios in other nuclei or mechanisms. For 
instance, it has been noted in that left-right symmetric models with hadronic right-handed currents have better 
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chances to be distinguished from v exchange via multi-isotope half-Ufe data. Extra-dimensional models with Kaluza- 
Klein neutrino states may also provide phenomenologically distinctive patterns |8i]. Moreover, additional information 
from angular and energy electron distributions in the final state might help breaking degeneracies (see, e.g., 0). 
The analysis of further (potentially more promising) mechanisms and observables, and of their relative discrimination 
in the presence of theoretical variants, is left to future works. 

A final remark is in order. Besides the QRPA framework adopted in this work, other independent approaches have 
been (and are being) adopted in the calculation of NME's for different nuclei, especially in the reference case of light 
Majorana neutrino exchange. Recent notable contributions include, e.g., applications of the Large Scale Shell Model 
(LSSM) ;435, of the Interacting Boson Model (IBM) [H], of the Projected Hartree-Fock-Bogoliubov method (PHFB) 
[45| . and of the Energy Density Functional method (EDF) [4^; see the review in [13]. The comparison of results 
obtained with these and other methods (or their variants) shows that the current spread of independent, state-of- 
the-art NME calculations exceeds the uncertainties estimated within the QRPA approach alone [47|, suggesting that 
further efforts are needed to constrain various theoretical approximations and systematics. These efforts, which would 
be boosted by future experimental observations of Oj^^/3 decay, will certainly help to understand better the source(s) 
of the degeneracy discussed in this work and, hopefully, to reduce the NME uncertainties at values small enough to 
allow the discrimination of at least two mechanisms. 



V. SUMMARY AND IMPLICATIONS 

We have considered nuclear matrix elements for neutrinoless double beta decay in four nuclei (^^Ge, ^^Se, ^™Mo, 
-'^'^^Te), within the self-consistent renormalized QRPA approach, for four Oi^2/3 decay mechanisms, mediated by light 
Majorana neutrino exchange {v)^ heavy Majorana neutrino exchange {N)^ i?-parity breaking supersymmetry {fi), and 
squark- neutrino {q). QRPA uncertainties have been evaluated by changing the model space size, the nucleon-nucleon 
potential, and the axial coupling. We have found that, within current QRPA uncertainties, the four mechanisms appear 
to be largely degenerate, and could co-exist in phenomenologically indistiguishable linear combinations (regardless of 
experimental errors and of possible complex phases). The degeneracy analysis is based on simple algebraic (rather 
than statistical) tools, which could be extended to NME estimates for other nuclei or mechanisms. 

As far as the degeneracy in the (^^Ge, ^^Se, ^°°Mo, ^■^'^Te) nuclei is concerned, one could, e.g., interpret their Oz^/3/3 
decay rates in terms of an effective phenomenological parameter of the form 

m'l^p ~rnpi3 + /^mpp, (21) 

where m/j^ is the usual effective Majorana mass for light v exchange, while Am^g^ embeds possible extra contributions 
from the (iV, ^, q) mechanisms, with either constructive (Am/3^ > 0) or destructive (Am^^ < 0) interference with 
the V mechanism. Note that, in the hypothetical case of nearly exact destructive interference (to^^ — 0), the process 
of QvPjS decay might be strongly suppressed and even escape observation in the four nuclei, despite the presence of 
Majorana neutrinos with mp^ > 0. 

A nonstandard degree of freedom AmpjS could be relevant in comparing 0i^/3/3 results with other searches for absolute 
u masses, such as /3-decay (probing m| = X^/il^ehpTi^) and precision cosmology (probing S = mi + TO2 -I- ms), 
especially if conflicting data emerge. An example is suggested by current results. Roughly speaking, for nearly 
degenerate neutrino masses (toi,2,3 — rn,y) it is E ~ Sm^, with typical cosmological constraints (S < 0.6 eV [i^) 
placing the upper bound ^ 0.2 eV. Assuming only the standard v mechanism, one has mp/s ~ m,yf (/ € [0.38, 1]) 
which, if mfjjs ~ 0.3 eV is inferred from the claim in [4^ |4l|, translates into a lower bound ~ mpp/ f > 0.3 eV. 
This well-known tension [48*1 could be relaxed by interpreting the OiyfSP claim, e.g., in terms of coexisting standard 
and nonstandard mechanisms with m^^ < 0.2 eV and Ampp ^JO.l eV, respectively (but adding up to m'^p ~ 0.3 eV); 
the same adjustment should then apply, within QRPA uncertainties, to the associated Oz^/3/3 signals in *^Se, "'^™Mo, 
and i30Te. 

We emphasize that our results refer to a specific set of four nuclei and four mechanisms, described within state- 
of-the-art QRPA nuclear structure calculations. Non-degenerate scenarios may well emerge for other 01/(3/3 nuclei, 
mechanisms and observables, or if prospective NME estimates in QRPA (or in alternative nuclear models) can reach 
significantly higher accuracies. We plan to explore some of these possibilities in future works. 
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Note Added. The NME estimates used in this work for various mechanisms and QRPA variants (see Table H]) are 
also being reported in another recent preprint '4^, where they are analyzed within different scenarios for coexisting 
mechanisms (with detailed discussions of non-interfering cases (SOj and/or of complex ?7's), as well as with different 
aims and perspectives. The approaches used in our paper and in (49| are largely non-overlapping, apart from the 
description of common NME inputs. The paper [49j analyzes what one can learn about the importance of the different 
mechanisms for the neutrinoless double beta decay, in particular, if accurate enough data and matrix elements will be 
available. The present paper shows that the current accuracy of the matrix elements does not prevent degeneracies 
and needs to be improved. One will also need future data with smaller errors for 2i//3/3 (to better constrain gpp) and 
precise enough half lives for Oz//3/3 decays to discriminate between different mechanisms for the neutrinoless decay. 
Finally, concerning non- interfering cases [sO] , it should be noted that the associated equations would be of the form 
T,i{Mil)-^r]j = isf in the unknown r]j [4^. Numerically, we find that the {MlY matrices and sub-matrices have 
reciprocal condition numbers about twice as large as the Mf (sub)matrices. However, since fractional differences due 
to different QRPA variants also double when Mf — (Mf)'^, there is no net gain, and the rjj solutions remain unstable 
within current uncertainties. 
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